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Thermodynamic ground states of platinum metal nitrides 
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The thermodynamic stabilities of various phases of the nitrides of the platinum metal elements are systemati- 
cally studied using density functional theory. It is shown that for the nitrides of Rh, Pd, Ir and Pt two new crystal 
structures, in which the metal ions occupy simple tetragonal lattice sites, have lower formation enthalpies at am- 
bient conditions than any previously proposed structures. The region of stability with respect to those structures 
extends to 17 GPa for PtNa. Calculations show that the PtN2 simple tetragonal structures at this pressure are 
thermodynamically stable also with respect to phase separation. The fact that the local density and generalized 
gradient approximations predict different values of the absolute formation enthalpies as well different relative 
stabilities between simple tetragonal and the pyrite or marcasite structures are further discussed. 
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Until recently none of the Platinum-Metal (PM) elements 
(Pt, Ir, Os, Ru, Rh, Pd) were known to form stable stable 
compounds with nitrogen. Several such compounds have now 
been synthesized primarily under conditions of high static 
pressure and temperature (in the range of 60 GPa and 2000K 
respectively) CL These include nitrides of Pt JlH, Ir flU, 
Os J3l and Pd lit] Except for PdN2 all these compounds have 
been shown to be at least metastable at ambient conditions. 
The resulting crystal structures have been investigated by sev- 
eral groups both experimentally and theoretically |2h8|], and 
by now consensus has been reached concerning the observed 
crystal structures and stoichiometry (one metal atom for every 
nitrogen dimer). PtN2 and PdN2 are formed in the pyrite crys- 
tal structure, a cubic phase with the metal atoms occupying fee 
sites. The nitrogen dimers are centered around the fee octa- 
hedral interstitial sites, oriented in all the four possible (111) 
directions, such that all the nearest neighbor dimers make an 
angle of 70.53° with each other. A rotation of dimers such 
that two pairs point in the [111] and [111] direction, respec- 
tively, results in the marcasite structure, the predicted ground 
state phase of RuN 2 , RhN 2 and OsN 2 Q. A small lattice 
distortion of marcasite then yields the monoclinic baddeleyite 
(CoSb 2 ) structure found in JjN 2 11 0] . 

The low-pressure phase diagrams of these systems are still 
not well understood. In particular their degree of thermody- 
namic stability has only been addressed briefly f^f], where it 
was proposed that the Pt-N compound is only metastable at 
low pressures. In this work, we address this issue by cal- 
culating the formation enthalpies of all the PM nitrides as a 
function of pressure. We also report the discovery of two new 
low-energy crystal structures that at low pressures are thermo- 
dynamically more stable than any of the crystal phases that 
have up to now been synthesized experimentally or calculated 
from theory. To our knowledge, these new crystalline phases 
have not previously been observed in any other compound. 

In spite of the apparent differences in the overall symme- 
tries of the pyrite, marcasite and baddeleyite structures, they 
all belong to the same class of crystal structures. They can 
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FIG. 1 : (Color online) Top: Projections of all possible distinct nitro- 
gen dimer orientations derived from the 12-atom conventional unit 
cell, along with the space group number and Hermann-Maugin sym- 
bol of the corresponding structure. Bottom: Energy/f.u. vs. volume 
curves from LDA for the above relaxed structures (solid lines) in the 
case of PdN2 . Each structure can be identified by matching the line 
color to the corresponding triangle color. Note that structure 7 corre- 
sponds to the black line. The green and red dashed lines correspond 
to the STaa and STab structures, respectively. 



be obtained by placing the PM atoms on an fee lattice and 
the center of (lll)-oriented nitrogen dimers on the octahedral 
interstitial sites of this lattice and allowing for atomic relax- 
ations. For each nitrogen dimer there are four possible orien- 
tations, and thus for crystals with the periodicity of the fee unit 
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cell having four distinct PM-N2 units, there are nine different 
crystal structures, two of which are the pyrite and marcasite 
phases. Note that as mentioned above, the baddeleyite struc- 
ture is merely a distorted marcasite structure. 

Geometrically each individual crystal structure in this class 
can be distinguished by the projections of its dimers centered 
at 0.5ax, 0.5ay, 0.5az and 0.5a (x + y + z) onto a (111)- 
plane. The dimer projections and space groups of the result- 
ing crystal phases are displayed in the top panel of Fig. Q] 
For reference, structures (3) and (8) correspond to marcasite 
and pyrite, respectively. To assess their theoretical equation- 
of-state we have performed density functional theory calcula- 
tions using projector augmented- waves as implemented in the 
Vienna ab-initio simulation package J9fl- Both atomic position 
and cell shape relaxations were performed. In fact, only struc- 
tures (8) and (9) maintained cubic symmetry after relaxation, 
while the others showed deviations of various magnitude. The 
resulting LDA energy vs volume curves for PdN2 are shown 
as the solid lines in the lower panel of Fig. [T] As expected, 
amongst these structures, pyrite and marcasite have the low- 
est energies in the low-pressure region. 

Besides the nine structures mentioned above, we also show 
in Fig. [T|the corresponding curves for two additional struc- 
tures (dashed lines). These are obtained by a continuous trans- 
formation from structure (2) by translating every other (001) 
plane in the [100] direction (corresponding to a zone-boundary 
phonon mode). By allowing the lattice to relax by a tetrago- 
nal distortion in the [001] direction, the metal ions become 
arranged in a simple tetragonal (ST) lattice. During this trans- 
formation the nitrogen dimers that initially were oriented in 
the (111) directions collapse into the (001) planes and now 
point in the (100) and (001) directions. In this way, a new 
class of crystal structures is created where metal nitride layers 
are stacked in different sequences. Each layer is composed 
of metal ions occupying square lattice sites and face-centered 
nitrogen dimers pointing along the edges of the squares and al- 
ways being perpendicular to the four nearest-neighbor dimers. 
Fig. |2] illustrates two particularly high-symmetry structures in 
this class with AA and AB stacking that have lower ener- 
gies than any previously proposed crystal structure for most 
PM nitrides at ambient conditions. The first structure, here- 
after denoted STaa, is simple tetragonal and belongs to space 
group 127 (P4/mbm). The metal ions occupy the 2a Wyckoff 
positions and the nitrogen atoms are located at the 4g posi- 
tions. The second structure, STab, is monoclinic and belongs 
to space group 12 (P2i /m). The metal ions occupy the 4e sites 
and the nitrogen atoms are located at the 4h and 4i positions. 

As is clearly seen in Fig. [U the ST structures at low pres- 
sures are significantly lower in energy than both the pyrite 
and marcasite structures. Also, their equilibrium volumes are 
shifted towards larger volumes and and the bulk moduli are 
smaller. Both LDA and GGA calculations yield small but fi- 
nite densities-of-states at the Fermi level for the ST structures 
1 1 111 . However, it is well-known that these approximations of- 
ten underestimate band gaps and e.g. predict Ge to be a metal. 
Hence, it is quite possible that the ST phases can be insulators 



3 



FIG. 2: (Color online) Conventional unit cells of STaa (left) and 
STab (right) viewed along the c-axis. Red, blue and yellow spheres 
represent metal ions, nitrogens in layer A and nitrogens in layer B, 
respectively. 



with small band gaps. 

The high-pressure PM-nitride synthesis results from the 
past few years have established the extent of (meta-)stability 
of many of these compounds. Since they have not been ob- 
served in nature, it is reasonable to rule out that they are ther- 
modynamically stable with respect to phase separation at am- 
bient pressure and temperature conditions. Based on exper- 
iments it is known that they are certainly stable at pressures 
above 50 GPa, and that PtN 2 (pyrite), IrN 2 (baddeleyite) and 
OsN 2 (marcasite), are at least metastable at zero pressure. 

There is however no real understanding either of the rela- 
tive stabilities among the compounds or of their their degree 
of metastability at low-pressure conditions. These are very 
important issues for the synthesis and potential application of 
these materials. As it is yet difficult to assess their equilibrium 
phase diagrams experimentally we have performed extensive 
first principles calculations of the thermodynamic stabilities 
of the six PMN2 compounds as a function of pressure. 

The thermodynamic stability of any PMN2 compound with 
respect to the separate phases as a function of pressure, is 
quantified in terms of the formation enthalpy: 



AH = Hp 



— HpM — Hn 



Under experimental conditions, nitrogen is molecular in 
its solid phase. At ambient pressure, it crystallizes into a 
hexagonal phase only at low temperatures and exhibits a poly- 
morphic phase diagram undergoing several structural trans- 
formations upon increasing pressure. It is thus quite difficult 
to perform accurate free energy calculations for this system. 
Therefore, the most reliable way of estimating the enthalpy of 
this system as a function of pressure is to split it into con- 
tributions from the dimer energy and the cohesive energy: 
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AH^. We note that the dominant contri- 



bution to i?N 2 at ambient conditions is due to the bonding en- 
ergy of the dimer, which can easily be calculated from first 
principles. We estimate the enthalpy of cohesion by calculat- 
ing it for the hexagonal phase at zero temperature and pressure 
within LDA/GGA and extract the change with pressure using 
the experimental equation-of-state data 112 1. 

Fig.[3]shows the calculated formation enthalpies within the 
LDA. For clarity, we consider here only the most important 
phases; pyrite, baddeleyite, marcasite and the two ST struc- 
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FIG. 3: LDA Formation enthalpies per formula unit of a) RuN2, b) RI1N2, c) PcOSb, d) OsN2, e) IrN2 and f) PtN2 as a function of pressure for 
marcasite or baddeleyite (MB), pyrite and the two simple tetragonal structures. 



tures. As baddeleyite and marcasite always display similar 
energies, we have chosen to only include the structure with 
the lowest formation enthalpy. 

We see that for the compounds formed with the group IX 
and X platinum metals, the ST structures become energeti- 
cally favored at lower pressures, while the group VIII nitrides 
form in the marcasite/baddeleyite (MB) structure at all pres- 
sures. The range of the stability of the ST structures versus 
pyrite and MB is consistently pushed towards lower pressures 
for the group IX nitrides as compared with the group X ni- 
trides. 

According to the LDA, the ST structures for PdN 2 and PtN 2 
are stable up to 13 GPa, while for IrN 2 and RhN 2 they are only 
stable at very low pressures. The high pressure crystal struc- 
ture of the group X nitrides is that of pyrite, while for all the 
earlier PM nitrides, the MB structures are the only structures 
observed at higher pressures. 

The same trend is observed in the GGA-calculations, dis- 
played in Fig. |4] except that the relative stabilities of the ST 
structures are more pronounced as they constitute the low- 
enthalpy phases up to 20 GPa for PdN 2 and PtN 2 and up to 
10GPaforRhN 2 andIrN 2 . 

A consistent feature observed in the phase diagrams is that 
the marcasite or baddeleyite structures are lowered in energy 
for the group VIII PM nitrides and the pyrite and ST structures 
are lowered in energy for the group IX and X PM nitrides. 

The formation enthalpies reported in Fig.[3]also contain in- 
formation on the absolute thermodynamic stability of these 
compounds. According to the LDA all but OsN 2 exhibit at 
least one stable phase with respect to the separate phases at 
ambient conditions. In fact, the lowest enthalpy phases for 



all considered PM-nitrides but OsN 2 are predicted to be sta- 
ble at all pressures. This is in contrast to the results from the 
GGA, which predict that the compounds are stable at 12 GPa 
for RuN 2 , 20 GPa for OsN 2 , around 15 GPa for RhN 2 , IrN 2 
and PtN 2 and 25 GPa for PdN 2 . 

In summary, we find from Fig.[3]the remarkable result that 
according to the LDA two novel crystal structures, never ob- 
served before, constitute the ground state phases of the group 
IX and X PM nitrides. PtN 2 is the most stable among all these 
compounds, with a cohesive energy of about 0.3 eV/f.u. In 
contrast, GGA predicts all the PM nitrides to phase separate 
at ambient conditions. For example, the PtN 2 ST phases are 
unbound within GGA by 0.6 eV/f.u. at zero pressure. They 
become thermodynamically stable at about 15 GPa. Upon fur- 
ther investigation of Figs. 0and[4]we fj nc j m at the difference 
between GGA and LDA enthalpies mainly stems from a rigid 
upward shift of the GGA curves by about 1.2 eV/f.u. This is 
not a surprising result as it is well known that LDA usually 
predicts larger cohesive energies and smaller lattice constants 
for solids than GGA. To compare with other nitrides, we find 
in the literature LDA formation energies of zinc-blende GaN, 
InN, A1N compounds that are 0.4 - 0.6 eV/f.u. lower than the 
corresponding values from GGA 1 1311 . This is less than half 
the difference found for the PM nitrides in this work. We at- 
tribute this to the presence of the nitrogen dimers. To prove 
this point, we calculated the formation energy of PtN in the 
zinc-blende structure and found that the LDA prediction is 
smaller than the GGA by about 0.4 eV/f.u. 

Hence, the difference between the LDA and GGA cohesive 
energies of the PM nitrides is consistent with previous litera- 
ture. As discussed above, this difference has important ram- 
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FIG. 4: GGA Formation enthalpies per formula unit of a) RuN2, b) RI1N2, c) PdN2, d) OsN2, e) LN2 and f) PflSb as a function of pressure for 
marcasite or baddeleyite (MB), pyrite and the two simple tetragonal structures. 



ifications for the low-pressure phase stability of these com- 
pounds. Generally, GGA is expected to provide better co- 
hesive energies and equilibrium volumes than LDA. This is 
true for the zinc -blende nitrides (GaN, A1N, InN) mentioned 
above. However, in the case of PM nitrides, LDA seems to 
slightly outperform GGA with regards to lattice constants. For 
example, LDA predicts the lattice constant of PtN2 to be 4.80 
A, in perfect agreement with experiment, as compared with 
the GGA lattice constant of 4.88 A 0]. Nevertheless, for the 
cohesive energies we believe that the GGA provides the most 
accurate values. 

When comparing with experiments we find that the calcu- 
lated formation pressures are less than half the observed syn- 
thesis pressures. For example, the measured synthesis pres- 
sure for the nitrides of Pt and Ir is approximately 50 GPa 
JH,^, and 58 GPa for palladium nitride Qj, while the calcu- 
lated formation pressures within GGA for the same systems 
are 17, 15 and 27 GPa respectively. At the experimental syn- 
thesis pressures, GGA predicts formation enthalpies of -1.4, 
-1.1 and -1.4 eV/f.u., respectively for these systems. We at- 
tribute this discrepancy to the presence of large kinetic barri- 
ers to the formation of the PM nitrides when synthesized in the 
diamond anvil cell. The latter may be associated with for ex- 
ample a strain due to lattice mismatch between the compound 
and the parent metal on which it is grown. The existence of 
large energy barriers to the formation of nitrides is well known 
in the synthesis literature. For example, according to GGA, 
GaN has a formation enthalpy of -1.7 eV/f.u. at zero pressure 
1 1311 but it requires a synthesis pressure of 2 GPa IU4I1 . 

In conclusion, the LDA clearly overestimates the stability 
of the PM nitrides. The GGA also predicts the transition pres- 



sures to be far below the experimental values. These results 
are however consistent with the large kinetic barriers found for 
zinc -blende nitrides. The only compound that exhibit a ther- 
modynamically stable ST-phase is PtN 2 . The combination of 
large barriers to formation and small region of stability make 
the synthesis of these phases in the diamond anvil cell diffi- 
cult. However, as we indicated earlier, these structures can be 
obtained through martensitic transformations from marcasite 
or pyrite structures. Hence careful experimental work at pres- 
sures in the window of stability of of the ST structures may 
be required to produce these phases. Furthermore, since the 
associated Raman spectra are probably weak another primary 
in-situ diagnostic may be required such as x-ray diffraction. 

This work was performed under the auspices of the U.S. 
Department of Energy by Lawrence Livermore National Lab- 
oratory under Contract DE-AC52-07NA27344. 
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